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Understanding of Ultra-Cold Neutrons Production in Solid Deuterium. 
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Our recent neutron scattering measurements of phonons and other quasi-particle excitations in 
solid deuterium (SD2) and the extraction of the density of states for phonons and rotational tran- 
sitions in SD2 have led us to a new understanding of the production of ultra-cold neutrons (UCN) 
in SD2. The UCN production rate reaches a maximum at an equivalent neutron temperature of 
T n =40 K for a neutron flux with Maxwellian energy distribution. The cross section for UCN pro- 
duction in SD2 has been determined by using the density of states Gi(E) in combination with the 
incoherent approximation as well as by a direct calibration of our measured neutron cross sections 
with the known cross section of the J = 1 1— > rotational transition in deuterium. Using this 
cross section we deduced the production rate of UCN in SD2 which agrees quite well with direct 
measurements of this energy averaged UCN production cross section. 
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Ultra-cold neutrons (UCN) are slow enough to be con- 
fined [l[ in traps, which can be formed by materials with 
a high Fermi potential (up to 300 neV) or by a magnetic 
field (60 neV/T). UCN can be observed for several tens 
of minutes in these traps, and are excellent tools for high 
precision measurements, as the life time of the neutron 
itself [2, 3], and the search for a possible electric dipole 
moment of the neutron [ij (current upper limit 2.9T0~ 26 
e- cm). 

Powerful UCN sources are needed for the here men- 
tioned experiments, and different groups [5-8] are work- 
ing on the development of strong sources, based on sD 2 
as an converter for down-scattering of thermal or sub- 
thermal neutrons into the UCN energy region (typical 
E <300 neV). A converter based on sD 2 should be op- 
erated at temperatures below T<10 K in order to avoid 
subsequent upscattering of UCN by phonons within solid 
deuterium [9]. The understanding of the downscatter- 
ing of neutrons into the UCN region crucially relies on 
a detailed knowledge of the energy loss of thermal or 
subthcrmal neutrons in the sD 2 converter material. One 
major reaction channel is based on phonons in the solid 
deuterium crystal, which are excited by the neutrons. 

We have recently measured the phonon system in 
sD 2 by neutron time-of-flight measurements at the IN4 
(ILL Grenoble) and at the TOFTOF (FRM II, Munich). 
These measurements are described in detail in another 
paper 0. 

Fig. [1] shows two examples of our neutron scatter- 
ing data (angular averaged) for two different concentra- 
tions c D of ortho-D 2 molecules. On the side of energy 



gain of the neutrons (positive E) we can identify the ro- 
tational transition J = 1 n- 0. The intensity of this 
peak decreases with reduced para-D 2 molecule concen- 
tration c p — 1 — c G . A similar rotational transition (the 
J = H > 1 transition) can also be observed on the neu- 
tron energy loss side. Furthermore it is remarkable that 
the elastic cross section (E =0 meV) of sD 2 is quite large. 
Our data indicates that the elastic cross section makes 
up approximately 50% of the total cross section at ther- 
mal neutron energy. Another interesting feature is the 
neutron scattering intensity on the neutron energy gain 
side close to the elastic peak (£"=0 meV). An increase 
of the concentration of the para D 2 molecules leads to 
a larger neutron up-scattering close to the elastic line. 
The origin of this cross section is very likely induced 
by phonons, which are interlinked with the rotational 
transitions J = 1 1-4 of the para molecules. This up- 
scattering has serious implications to the up-scattering 
of UCN in the sD 2 converter material, and reduces the 
achievable density of UCN in sD 2 . A detailed analysis of 
these effects will presented in a forthcoming paper. 

In addition the scattering function S(Q,E) is shown 
for fixed values of momentum transfer Q in fig. [2] 
The dynamical scattering function S(Q,E) comprises 
the phonon branches of the hcp-sD 2 [11[ crystals and 
the J = M> 1 rotational transition. The acoustical 
phonons (transversal and longitudinal) are clearly visible 
in S(Q, E) while the optical phonon branches are not. 
The rotational transition J = M> 1 shows a significant 
Q - dependence and it is purely incoherent. The posi- 
tion the J = 0h>1 transition peak (Eq\=7A meV) does 
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FIG. 1: An example of a dynamical neutron cross section of 
solid D2 at T — 7 K. Comparison of two ortho concentra- 
tions Co =66.7% (□) and c D =98% (CO- Data from TOFTOF 
measurements at the FRM II. Initial energy of the thermal 
neutrons is Eo=20.4 meV. 



neutron cross section -j^ (averaged over the scattering 
angle, thus Q). Vice versa it is also possible to extract 

using Turchin's theory 
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G\(E) from a measured 

[l2| for this cross section applying the incoherent approx- 
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imation. The absolute normalization of 

needed, because the extracted G\(E) has to be normal- 
ized to unity anyway. This approach is only valid within 
the Born-approximation. 

The method of determination of G\ (E) from our data 
is described in [lfj in detail. Contributions of higher 
order multiphonons to -j^j are incorporated. 

The result of our analysis [l(| concerning G\(E) and 
a comparison with a Debye model is shown in fig. |4] 
The characteristics of G\(E) was already discussed in 
detail in [l(| but it is worth to summarizing these results 
again, because they have an important impact on the 
UCN production cross section. One major feature is 
the occurrence of these excitations above E^IO meV in 
the one-quasi-particle density of states G\(E). These 
excitations above E~10 meV are depending on the 
concentration of ortho-D2. The optical phonons in the 
region of E [8-10 meV] are not clear visible in G\(E). 



not change much at different Q- values. The cross section 
d(Tj=o^i/dQ for this transition is shown in fig. [3] This 
cross section was extracted from our data by a integration 
of S(Q,E) at £ni=7.4 meV with a width of AE= lmeV 
at different Q-values. The Q-dependence of this cross 
section was fitted using the incoherent scattering length 
of the J = i-> 1 transition jig ] 



J=0<-^1 



(Q)~fi(Qa a /2) 



- L 3 Q 2 (u 2 ) 



(1) 



The parameter a s = 0.74 A is the distance of the 
deuterons within the D2 molecule, while (u 2 ) is the mean 
square displacement of the D2 molecule in the lattice. 
The form factor jx for the J = t->- 1 transition is de- 
scribed by a spherical Bessel function of first order. The 
result for (u 2 ) ~ 0.245 ± 0.02 A is in good agreement 
with earlier published results ll|. The peak-like devia- 
tions in fig.|3]at certain Q-values are the result of coherent 
neutron scattering on phonons. 

With the aid of our neutron data we are able to 
determine the UCN production cross section by two 
ways. One way is the determination of the quasi-particle 
(phonons and rotational excitations of the D2 molecules) 
density of states G\(E). The other way is the direct 
integration of the dynamical neutron cross section in 
the kinematical region along the free neutron dispersion 
parabola. 

With the knowledge of the quasi-particle density of 
states G\(E) it is possible to calculate the dynamical 



In the case of UCN production the energy transfer of 
the down-scattered neutron E — Eq — Ef is approxi- 
mately equal to the initial neutron energy Eq (Ef = 
E v « E ,Eu: UCN energy). 

The total cross section for UCN production can be 
calculated by 



o'ucn(Eq) 



da(E ] 
dE 



dEr, 



(2) 



and is shown in fig. [5l It can be compared with recent 
published data of the UCN group at the PSI [l3| which 
has studied UCN production at a cold neutron beam 
(Eo ~1.4 meV to 20 meV). The agreement of data and 
our calculation, using our data for G\(E) and Turchin's 
incoherent approximation is surprisingly good. The cal- 
culated cross section comprises the contribution of one- 
quasi-particle and two-quasi-particle excitations. Three- 
quasi-particle excitations do not appear below E~14 meV 
(see fig. 9 in [13) and are rather small in their contri- 
bution. The authors of 13[ compared their results to a 



calculated cross section based on a multi-phonon Debye 
model. They also reported, that using a 'more realistic 
model' for the density of states ( Yu et al. 11, TtJ) leads 
to a considerably worse fit of their data. Using a Debye 
model with a cut-off at Ub-Qd =9.5 meV (Debye temper- 
ature Od=110 K) only the contribution of two-phonon 
excitation have to be included in order to fit the measured 
cross section of 10 -7 barn at E=14.7 meV 

13j |. This also applies for the density of states model 
of Yu et al. Our own calculated cross section shows a 
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FIG. 3: Cross section da j— 0^1/ dil of SD2 for c =95.2% 0-D2. 

■ - Data from IN4 measurement. Black line - Fit result using 
FIG. 2: S(Q, E) of sD 2 for c D = 95.2% at T = 4 K for different eq. [Q 
fixed values of Q. Data from IN4 measurements. 



combination of one-quasi-particle and two-quasi-particle 
excitations contributing at these energies. In detail our 
data indicate (see fig. [5]) that the UCN production cross 
section is mainly determined by one-quasi-particle exci- 
tation for energies below E=15 meV ( see fig. [4]). Even 
though the two-quasi-particle contribution can not be ne- 
glected in the region of E 5-25 meV. 

The application of the incoherent approximation in the 
case of sD 2 has certainly to be questioned. The sD 2 crys- 
tal scatters neutrons more coherently than incoherently. 
This fact leads to the question: Is it possible to get the 
UCN production cross section directly from the neutron 
scattering data? 

The easiest way of determining the cross section for 
UCN production is the use of the dynamical scatter- 
ing function S(Q,E = ^Q 2 ) in the (Q,E)-phase space 
along the free neutron parabola (dispersion curve with 
E ~ Eq = ^k" 2 and Q ~ k ). This approach equally 
rely on the Born approximation. The UCN production 
cross section can be determined by 



&ucn(Eo) 



Y S{ko ^ ko) 3 ku Eu ■ 



(3) 



The term E ~ En = ^—k 2 is the energy of the incom- 
ing neutron in the down-scattering process, while k™ ax 
and Ey ax are the upper limits for the UCN momentum 
and energy. In order to obtain absolute cross sections 
the dynamical scattering function of solid deuterium ex- 
tracted from our neutron scattering data has to be cali- 
brated to absolute values. This can be done by using the 
rotational transition J — 1 M> as reference process for 
calibration. The cross section for the J = 1 n- transi- 



tion was calculated by Hamermesh and Schwinger 14 1 
for ortho- and para- deuterium molecules (gas). These 
calculation include also thermal movements (maxwellian 
distribution of velocities) of the deuterium molecules in 
the gas. The application of these calculations to the rota- 
tional transitions of the D 2 molecules, which are pinned 
in a crystal at fixed positions may of course be question- 
able. However the rotational transitions in sD 2 are not 
very much hindered by the crystal binding of the crys- 
tal [HI. Indeed the molecules are of course not at rest, 
they perform motions around their equilibrium positions. 
These movements are caused by the zero-point-motions 
and also by existing phonon excitations in the crystal. 
A similar problem was considered by Lamb [lrij in the 
past. He calculated the capture of neutrons by atoms in 
crystals. His calculations for crystals with weak lattice 
binding delivered the finding, that atoms in such crystals 
can be treated as a gas with an effective temperature, 
which depends on the Debye temperature of the crystal. 
This effective temperature can be higher than the real 
temperature of the solid and is just a result of the aver- 
age energy of the vibrational degrees of freedom plus the 
zero-point-motion. We used this ansatz to calculate the 
effective temperature (T e ff ~40 K at a deuterium tem- 
perature of 4 K) of the D 2 molecules and used it for the 
calculation of the rotational transition (dj—i^o) of our 
quasi free molecules, following the theory of Hamermesh 
and Schwinger. 

The cross section aj— 1^,0 i n pure para deuterium has 
a value of 07=1,^0=0.61 barn for E =17.2 meV ( corre- 
sponding to the energy of the incoming neutrons at the 
IN4 experiment) at a deuterium temperature of T=4 K. 
The comparison of our neutron scattering data of natu- 
ral deuterium (co=66.7%) and ortho enriched deuterium 



4 



0.32 



0.24 



> 



0.16 



0.08 



0.00 




2 4 6 8 10 12 14 16 18 20 
E [meV] 




10 15 
E [meV] 



FIG. 4: Comparison of the one-quasi-particle density of states 
of solid D 2 for c D = 66.7% (O) and c = 98% (□) at 
T = 7 K with a Debye model (line). Data from TOFTOF 
measurements (energy resolution: AE ~ 1.24 meV). The 
SD2 poly crystal is prepared by fast freezing (several min- 
utes) of liquid D2 down to a temperature of T <10 K. 
J °° GDOS(E) ■dE=l 



FIG. 5: UCN production cross section of c = 98% solid D2. 
UCN energy range 0-150 neV inside the solid D2. Solid line 
- cross section calculated in incoherent approximation, using 
Gi(E) from Dashed line: one-quasi-particle contribu- 

tion. Dotted line: two-quasi-particle contribution. □ - data 
from measurements at the PSI 



indicated a ortho concentration of co = (95.2 ± 0.3)%. 

Scaling with the para concentration of our sam- 
ple (cp = 1 — co =4.8%) thus leads to a value of 
cj = ih^o =0.029 barn for our ortho enriched solid D2 sam- 
ple. This value for uj = \^q can be used to calculate a 
scaling factor for S(Q, E) = k ■ S dat a(Q, E), 

By means of this value it is possible to calculate the 
total neutron scattering cross section for Eq=17.2 meV 
neutrons: atot{Eo=17 .2 meV)=23.4 barn. The assump- 
tion is of course only valid, if the kinematical area of our 
neutron scattering experiment covers most of the possi- 
ble scattering processes. A good cross check of this value 
is the calculated value for the total cross section for solid 
poly-crystalline deuterium, using the incoherent approx- 



imation of Turchin 12 



The result of this calculation leads to a value of 
(Ttot(Eo— 17 .2 meV)=23.8 barn, which is close to the 
value of our calibration. The dynamic response func- 
tion S(Q,E), resolved from our data analysis, was not 
corrected for multiple scattering effects. Therefore the 
estimated error of this calibration is approximately 12% 
[20l | . The result of this calibration and determination of 
the UCN production cross section as function of the the 
energy of the incoming neutrons, and a comparison with 
measurements of this cross section 13j is shown in fig. 
[6j The agreement of the absolute values of the measured 
cross section with the resolved values form our S(Q 7 E) 
data is reasonable. The comparison of the calculated 
UCN production cross sections, extracted from the inco- 



herent approximation and parabola method shows (see 
fig. [5] and fig. [5]) a discrepancy in the region of E ~ 
6 meV. The cross section, determined by the parabola 
method exhibits a pronounced maximum in the region of 
E ~ 6 meV as compared to the incoherent approximation 
result. 

In fig. [7] the dynamical scattering function S(Q,E) 
(neutron energy loss side) of fast frozen solid ortho- 
deuterium is shown. The black line corresponds to the 
dispersion of the free neutron. 

The parabola of the free neutron crosses the acousti- 
cal phonon dispersion curve at E ~ 6 meV(see fig. [7J. 
At this point, the UCN production cross section is pre- 
dominantly determined by coherent scattering. This can 
explain a deviation from the production cross section 
in incoherent approximation. Nevertheless the general 
agreement of the incoherent approximation with the PSI 
data is remarkable (see fig. [SJ 

The effective UCN production rate P (UCN cm -3 s" 1 ) 
from a neutron beam is determined by integrating the 
product of the UCN production cross section da/dE and 
the spectral flux d$/dE (Maxwell spectrum with effec- 
tive neutron temperature T n ) of the incoming neutrons 
over the finite UCN energy range (E™ ax maximum al- 
lowable UCN energy) and over the energy spectrum of 
the incoming neutron, still assuming E = Eq. 



P(T n ) = N D2 



d$(T n ) da 

-^-—tEvdE*. (4) 
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FIG. 6: 

UCN production cross section solid D2 of c = 95.2%. A 
UCN energy range of 0-150 neV inside the solid D2 is 
assumed. Cross section determined by an integration of 
S(Q, E) along the free dispersion of the neutron (sample - 
fast frozen solid deuterium (T =4 K)); data from IN4 
measurements. Blue □ - Eq =17.2 meV. Red filled O " 
So =67 meV, 118j. ■ - direct UCN production data from 
measurements at the PSI [l3| • The PSI data do not show a 
distinct peak at E ~ 6 meV 
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FIG. 7: 

S{Q,E) of sD 2 for c D = 95.2% at T = 4 K. Data from IN4 
measurements. Black parabola - dispersion of the free 
neutron. 



In fig. [5]the result for the UCN production rate in solid 
ortho deuterium, exposed to Maxwellian shaped neutron 
flux for different effective neutron temperatures is shown, 
using eq. HI 

The main conclusion from these results is the new 
understanding of possible higher energetic loss channels 
(one-quasi-particle and two-quasi-particle excitations) in 
solid deuterium for the down scattering of thermal or 
cold neutrons in the conversion process to UCN. The 
best value for the effective neutron temperature is in 
the region of T n —40 K, which is larger than the value 



FIG. 8: Calculated UCN production rate of c a = 98% solid D 2 
for different Maxwellian neutron spectra with effective neu- 
tron temperature T„. UCN energy range - 0-150 neV inside 
the solid D2. Neutron capture flux $0 = 1 ■ 10 14 cm -2 s _1 . 
Solid line - total production rate (one- and two particle ex- 
citation). Dashed line - one-particle production rate. Doted 
line - two-particle production rate. 



(T„ ~30 K) reported by Yu et al. [T7| in an earlier pub- 
lication. A similar result was obtained by Serebrov et 
al. [HI using a Debye-model for SD2 in a theoretical cal- 
culation, based on the incoherent approximation. The 
agreement of our result with the results of the calcula- 
tion applying a Debye-model is not a surprise, because 
they include mulit-phonon excitations in their calcula- 
tions, while Yu et al. considered only one-phonon exci- 
tations. 

In summary, new neutron scattering data of solid 
deuterium leads to a better understanding of UCN 
production in this converter material. The new results 
for the density of states in SD2 and the results for the 
UCN production cross section, extracted directly from 
the dynamical scattering function S(Q, E) predict a 
significant UCN production cross section for incoming 
neutrons with energies higher than Eq > 10 meV. 
This observation was also confirmed 13[ by direct 
measurements of the UCN production cross section. An 
optimized sD 2 UCN source should be exposed to a cold 
neutron flux with an effective neutron temperature of 
T n -40 K. 
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